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Abstract: Grazing exclusion is one of the most efficient approaches to restore degraded grassland but 
may negatively affects the recovery of species diversity. Changes in plant species diversity should be a 
consequence of the ecological assembly process. Local community assembly is influenced by 
environmental filtering, biotic interactions, and dispersal. However, how these factors potentially 
contribute to changes to species diversity is poorly understood, especially in harsh environments. In this 
study, two management sites within a S/pa breviflora desert steppe community (typical natural steppe) were 
selected in northern China. In one of the two management sites, grazing has been excluded since 2010 
and in the other with open grazing by sheep. In August 2016, three plots were established and 100 
sampling units were created within each plot in a 5 mX5 m area at the two management sites. To assess the 
effects of grazing exclusion on S. breviflora steppe, we analyzed the vegetation biomass, species diversity, 
soil organic carbon, and soil particle size distribution using paired T-tests. In addition, variation 
partitioning was applied to determine the relative importance of environmental filtering and dispersal 
limitation. Null mode analysis was used to quantify the influence of biotic interactions in conjunction with 
EcoSim niche overlap and co-occurrence values. Our results demonstrated that (1) species diversity 
significantly decreased and the main improvements in soil quality occurred in the topsoil 0—10 cm after 
the grazing exclusion; (2) environmental filtering was important for community assembly between grazed 
and fenced grassland and this appears particularly true for soil particle size distribution, which may be well 
correlated with soil hydrological processes; and (3) however, competitive exclusion may play a significant 
role within the exclusion. The multiple pathways of assembly may collectively determine negative effects 
on the restoration of species diversity. Therefore, designers should be aware of the risk of reducing 
grazing exclusion-induced species diversity and account for manipulating processes. This in turn will 
reduce dominant species and promote environmental heterogeneity to maximize species diversity in 
semi-arid regions. 
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1 Introduction 


Grasslands have experienced dramatic degradation under human disturbances and global climate 
change, both of which threaten biodiversity and ecosystem functioning (Bai et al., 2012; Shao et 
al., 2016). Recently, many ecological engineering approaches have been used to restore degraded 
grassland such as shrub planting, grazing exclusion, species introduction, conversion of land use, 
or altered grassland management. Biodiversity is critical for ecosystem resilience 
(Coiffait-Gombault et al., 2012). Thus, the ultimate goal of ecological restoration should be to 
maintain biological diversity and the associated ecological services (Cadotte et al., 2011; Peco et 
al., 2012; Wu et al., 2014). Species diversity is a component of biodiversity that should be highly 
correlated to ecosystem properties such as resource use, primary production and stability 
(Balvanera et al., 2006). Therefore, restoration of plant species diversity is crucial for degraded 
grassland ecosystem and has attracted significant attentions in recent years. Grazing exclusion is a 
simple way to induce self-recovery and has been regarded worldwide as one of the most efficient 
approaches to achieve this (Wu et al., 2009; Wu et al., 2010; Jing et al., 2014). However, grazing 
exclusion may have negative effects on plant diversity (Jing et al., 2014). Although the negative 
effects of grazing exclusion on species diversity are well documented, the processes that underlie 
variations in species diversity and how this relates to natural vegetation restoration are not well 
understood. 

Species diversity, measured as species richness, abundance, and Shannon diversity index, is 
regarded as an outcome of species interactions or community adaptation to the biotic or abiotic 
environment in conjunction with evolutionary factors (Rice and Westoby, 1983). The alteration of 
species diversity is a consequence of the ecological assembly process (Mayfield et al., 2010). The 
tule of ecological assembly describes community structure and composition by passing them 
through the following ecological filters: the external abiotic environment (e.g., topography, soil, 
and precipitation), internal dynamic including biotic interactions, and dispersal (G6tzenberger et 
al., 2012; Nathan et al., 2016). Indeed, ecological assembly theories, advances, and challenges 
have been extensively reviewed (Silvertown, 2004; Vellend, 2010; Weiher et al., 2011; 
Gotzenberger et al., 2012; Hillerislambers et al., 2012), and the mechanisms of community 
assembly that are influenced by processes operating in diverse, spatiotemporal ranges are 
well-known. Species found over continental scales are constrained by evolutionary history and 
biophysical conditions (Virtanen et al., 2006); whereas, those found within local scale community 
assemblies reflect external variables and internal dynamics (Irvine et al., 2009). The intra and 
interspecific interactions within small spatial scales that have only a few centimeters or 
millimeters between neighbors are likely to be influenced by biotic filters (G6tzenberger et al., 
2012). The community assembly processes are also disturbance dependent. Grazing is the most 
common disturbance regime that has a pronounced effect on species assembly (De Bello et al., 
2013). Grazing directly influences competition and introduces heterogeneity into grasslands via 
selective grazing, treading, seed distribution, and nutrient redistribution (Jerrentrup et al., 2015). 
For example, selective foraging causes a local decrease in the competitive exclusion of species, 
thereby allowing subordinate species with low competitive abilities to coexist and increase 
species richness (Grime and Mackey, 2002; Mayfield et al., 2010). Moreover, trampling by 
grazing animals and subsequent fine-scale soil disturbances may help to create niches for 
gap-colonizing species (Hofmann and Isselstein, 2004). In addition, patches of bare soil lacking 
vegetation cover may offer space for the propagation and dispersal of seeds (Jerrentrup et al., 
2015). In contrast, grazing exclusion can be followed by a spontaneous recovery of soil processes, 
disturbance regime, and plant-plant interactions (Coiffait-Gombault et al., 2012; Wu et al., 2014; 
Wang et al., 2015). Shifts in species coexistence patterns and community assembly processes may 
occur after grazing pressures are released through the use of mesh fencing. For example, 
long-term grazing exclusion enhanced the litter layer accumulation, competition for water and 
light resources, and competitive exclusion between different species (Guo, 2005; Wu et al., 2009). 
In contrast, a lack of grazing in grasslands may improve niche differentiation (Mason et al., 
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2011), which promotes species coexistence and community diversity (Siepielski and Mcpeek, 
2010; Thompson et al. 2010; Freschet et al., 2011). In addition, internal self-organization is a key 
driver for community assembly; whereas, external variables play a less significant role in 
community assembly after grazing pressures are released over a long term (Lawley et al., 2013). 

Natural grasslands cover 40% of area in China (Hua and Squires, 2015), making it the largest 
ecological system in the country with important ecological functions for ecological security, 
including water and soil conservation, biodiversity, wind resistance, sand fixation, and air 
purification (Shao et al., 2016). Approximately 80% of the grasslands in the regions of western 
China have undergone the worst degradation in the world as a consequence of overgrazing and 
other human activities (Wu et al., 2014). Restoration, followed by community self-organization in 
arid and semi-arid grassland ecosystems, has complex outcomes driven by internal interactions 
and does not follow a linear, predictable trajectory (Lawley et al., 2013). Importantly, the 
trajectories of community development may not necessarily lead to the desired state (Andrade et 
al., 2015). Previous studies that focused on grassland restoration have examined changes in 
vegetation (e.g., composition and structure, diversity, and biomass) and soil characteristics (e.g., 
water content and carbon storage) in different grassland type, primarily focus on restoration 
evaluation (Zhang et al., 2011; Jing et al., 2014). These studies reported that grazing exclusion in 
semi-arid grasslands significantly improved vegetation with higher species diversity in the early 
stages of restoration, but decreased species diversity and biomass with increasing time of 
restoration. This may be due to the long-term grazing inducing litter layer accumulation, 
competition for water and light resources, and competitive exclusion between different species 
(Wu et al., 2009). Recent studies have focused on patterns of beta diversity and alteration in 
ecosystem functioning in China's grasslands (Bai et al., 2012; Chi et al., 2014). However, little is 
known about the relative contributions of environmental filtering, biotic interactions and dispersal 
on community assembly at a local scale. For the maintenance and restoration of diversity, 
understanding the relative importance of different factors in the process of community assembly 
at a local scale is necessary and required for establishing guidelines for restoration of degraded 
grassland eco-environments (Lawley et al., 2013). In this study, two management sites of Stipa 
breviflora steppe (typical natural steppe) were selected in a desert steppe ecosystem in northern 
China. One was a fenced grassland where grazing had been excluded since 2010, and the other 
was a sheep-grazed grassland at a level of grazing intensity of 1.5 sheep/hm?). Because S. 
breviflora is a typical species in continental steppes, it is regarded as a mature ecosystem that can 
be considered a reference system for restoration (Bai et al., 2012; Jing et al., 2013, 2014). The 
objectives of this study were (1) to investigate the effect of grazing exclusion on vegetation 
characteristics and soil quality; (2) to assess the relative contributions of environmental filtering, 
biotic interactions, and dispersal limitation on community assembly between the two management 
sites; and (3) to compare the shifts in the different processes in order to gain insight into 
appropriate strategies for degraded grassland restoration. 


2 Materials and methods 


2.1 Study area 


The experiment was conducted in Maerzhuang village (37°27'30"—37°37'30"N, 106°37'30"— 
106°56'15"E; 1450 m a.s.1.) in Yanchi County, Ningxia Hui Autonomous Region, China. The area 
is located in the southwestern fringe of Mu Us Sandy Land and in the mid-western regions of the 
Loess Plateau. Harsh climatic conditions in this area make up the majority of the underlying 
forces causing land desertification. However, overgrazing is also known to exacerbate 
desertification. The annual mean temperature is 8.1°C, ranging from —8.7°C in January to 22.4°C 
in July. The mean annual precipitation is 295.1 mm, mainly occurring in June to September (Tang 
et al., 2015). The mean annual potential evaporation is 2136 mm. The main soil types in this area 
are sierozem, loess, and orthi-sandic entisols, all of which have low fertility and loose structures 
(Liu et al., 2014). The representative vegetation include S. breviflora, Leymus chinensis, and 
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Agropyron cristatum communities. However, the representative vegetative communities have 
been greatly altered, primarily due to the long-term overgrazing. 


2.2 Experimental design and data collection 


The management status of the studied grasslands in the region ranges from fenced grasslands (Fig. 
1), where grazing has been excluded since 2010, to grasslands that are grazed by domestic 
animals (mainly sheep). In late July 2016, we selected two study sites for sampling: (1) inside the 
fenced grassland, where livestock had been excluded for six years (from 2010 to 2016); and (2) 
outside the fenced grassland, where grazing was allowed all year round. The predominant plant 
species at the two sites were S. breviflora (perennial grass), and Cleistogenes squarrosa. Three 50 
mx50 m plots (500 m apart from each other) were selected at both the grazed and fenced sites. A 
total of six plots were established and each plot was divided into three parts. 


Inside the fence 


The fence line 


Fi 
he 
7 


Outside the fence 


Fig. 1 The landscape of study sites. The middle picture was taken using unmanned aerial vehicle 
(DJIS1000+Canon5D Mark III). The fence line is clearly visible with a sharp contrast in the vegetation density 
and cover on its both sides. 


2.2.1 Sampling for species-area relations 


There were 13 subplots (0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 0.64, 1.28, 2.56, 5.12, 10.24, 20.48, 
and 40.96 m°?) created within a nested quadrat for measurements of species-area curves. Species 
presence and absence were recorded in all 13 subplots. This provided information in determining 
the minimum sampling scale in Section 2.2.3 (Fig. 2). 

2.2.2 Sampling for vegetation and soil inside and outside the fence 

Five randomly distributed subplots (1 mx1 m each) were established within each 50 mx50 m plot 
to measure vegetation and soil properties. Above-ground biomass of green plants within each 
subplot were harvested, and the number of species, plant height, cover, abundance, and frequency 
for each plant species were recorded. The specific operations for a detailed vegetation survey and 
biomass measurement were given by a study in the same region (Tang et al., 2015). At the same 
time, three random soil samples in each quadrat were collected with a 5-cm diameter soil auger 
from four depths: 0-5, 5—10, 10—20 and 20—40 cm. Soil samples at the same depth in each subplot 
were then mixed to a single sample. 
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2.2.3 Sampling for variation partitioning and the pattern of species coexisting inside and outside 
the fence 


One 5 mx5 m subplot was established in each 50 mx50 m plot for assessing variation in 
community assembly between the grasslands inside and outside the fence. Vegetation surveys 
were carried out within the subplots, which were further divided into 100 sampling units (50 
cmx50 cm). Trampling from domestic animals can cause soil compaction, which alters soil 
structure, soil bulk density, and the transport of oxygen and water. These changes may be 
detrimental to plants in grasslands (Van Klink et al., 2015). Soil compaction and water content 
can be detected directly in the field using a stratameter and time domain reflectometry, 
respectively at a depth of 0-10 cm. At the same time, soil samples were taken at 0—10 cm depth. 
We collected a total of 100 soil samples and established 100 vegetation sampling units within 
each 5 mx5 m subplot. The plot scales, 0.5 and 5.0 m°, were used for assessing variation in 
community composition because: (1) the species-area curve saturates at 5.0 m? and the minimum 
size of sampling unit should be 0.5 m? (Fig. 2); and (2) co-existing species divergence may occur 
at the finest spatial scale. At those scales, individuals from different species are more competitive 
with neighbors for resources and space (Luzuriaga et al., 2012). For example, co-existing species 
at 50 cmx50 cm showed a convergence pattern caused by environmental filtering in dry 
semi-natural grasslands (De Bello et al., 2013). 
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Fig. 2 Species-area curves inside and outside the fence. Two saturation points at plot size 0.5 and 5.0 m? were 
determined, respectively. 


2.2.4 Laboratory analysis 

In the laboratory, soil samples were air-dried and passed through a 2-mm sieve. The roots and 
debris were carefully picked out and discarded. The samples were then divided into three parts. 
One part was stored at room temperature. The second part was passed through a 0.15-mm sieve 
for the analysis of soil organic carbon by dichromate. The third part was used for the particle size 
distribution, which was described in relative percentage of clay (<0.002 mm), silt (0.05—0.002 
mm), fine sand (0.1—0.05 mm), or coarse sand (2—0.1 mm). This analysis was determined using a 
Mastersizer 3000 particle size analyzer (Malvern Instruments, Malvern, England). 


2.3 Data analysis 


2.3.1 Species diversity, biomass and soil quality changes 

We used the average mean of the plots to determine the differences between the treatments (inside 
and outside the fence). Vegetation and soil properties were analyzed using paired-sample t-tests in 
SPSS software. Biomass measured by Grass, Leguminosae, Compositae and weeds, respectively. 
Plant diversity was assessed using the Richness index (R), Shannon-Wiener index (H’), Simpson 
index (D), and Pielou index (J). Plant diversity indexes were calculated using the "vegan" library 
of the R software (Oksanen et al., 2016). 

2.3.2 Relative roles of spatial factors and environmental filtering 


Based on the mean values of species diversity matrix for each site, we employed two methods to 
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disentangle the relative importance of the environmental and spatial process independently or in 
combination. The mantel and partial Mantel test, based on the environmental differences (as a 
proxy for environmental filtering) and spatial distance (as a proxy for dispersal limitation), were 
used to detect whether the species diversity depended on the environmental and spatial distance, 
which was based on Euclidean distance (Wang et al., 2013). 

Variation partitioning was then used to examine the relative contribution between 
environmental filtering and dispersal limitation based on Hellinger-transformed species data and 
standardized environmental data (Wang et al., 2013). First, the spatial eigenvectors (as a proxy for 
spatial factors) were obtained through the Principal Coordinates of Neighbor Matrices (PCNM) 
analysis across the coordinates of all samples. For a complete description of the method, see Dray 
et al. (2006). In a region where soil is nutrient deficient, some plants induce the build-up of 
organic matter and buffer against abiotic factors. Thus, their location may be more important in 
determining future spatial patterns than the soil substrate (Butterfield et al., 2010). At our study 
sites, the soil was of low fertility and the predominant species, S. breviflora had higher individual 
cover, density, and productivity in comparison with the other species. In this type of environment, 
S. breviflora is more likely to improve environmental conditions or enhance the survival and 
growth of the other species than in open spaces (Luzuriaga et al., 2012). Environmental variables 
tested included soil characteristics (organic carbon, water content and particle size) and 
aboveground biotic factors (cover and density of S. breviflora). Since there was strong 
collinearity between some environment variables, we first removed environment variables that 
were highly correlated and then used forward selection ("forward.sel" function in the R package 
version 0.0-8.0) to retain environment variables with a significant effect on plant diversity (Wang 
et al., 2013). Similarly, forward selection was also performed for the PCNM variables. The 
significant environmental and the PCNM variables (a=0.05, Monte Carlo permutation, n=999) 
retained in the final model were then applied to partition variation in plant diversity into four 
fractions: [a]=variation explained by the environmental factors independent of space, 
[b]=variation explained by the spatially structured environmental factors, [c]=variation explained 
by spatial factors independent of any environment factors, and [d]=residual variation (Oksanen et 
al., 2016). All analyses were conducted in R software. 


2.3.3 Pattern of species co-existing 


To determine if the communities inside and outside the fenced grasslands were randomly or 
non-randomly structured, we used null model systems to analyze species co-occurrence patterns 
based on the mean values of the presence-absence matrices for each site (McGee and Eaton, 2015; 
Miguel et al., 2016). Random patterns were compared with those observed in the community to 
obtain information about the underlying mechanisms (abiotic filters and biotic interactions) in 
community composition. To quantify patterns of species co-occurrence and identify the relative 
strength of these relationships, we applied a standardized effect size of the co-occurrence index in 
EcoSim 7.0 simulation software (Gotelli and Entsminger, 2001; Zhang et al., 2009). For a 
complete description, see Gotelli (2000). Assuming a normal distribution of standardized effect 
size (SES-C), 95% of the observation was located between —2.0 and 2.0. For example, for an 
SES-C score, an observed value of 1.5 or larger indicates nonrandom segregation and suggests a 
significantly less co-occurrence than expected by chance among species. This model 
demonstrates whether species co-occurrence in a community may be more competitive for space 
and resources; however, the opposite can be said when a SES-C score is less than —1.5 (Mcgee 
and Eaton, 2015). The P-values in the report of EcoSim 7.0 provided information about the 
Statistically significant difference of the observed index in relation to those obtained from 50,000 
simulations (Gotelli and Entsminger, 2001). We performed the fixed-fixed null model, which has 
a low risk of Type I error, to simulate the randomly constructed assemblages (Gotelli and 
Entsminger, 2001; De Miguel et al., 2016). In this null model, both row sums (species) and all 
columns (plots) remain constant and equal to the original data matrix. 

The niche overlap analysis was performed in EcoSim 7.0 to determine if the communities 
inside and outside the fenced grasslands were structured around an environmental condition or a 
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limited resource. The niche overlap model tests variation among species in a niche utilization. We 
selected the randomization algorithm conditions to exam the possibility that changes in certain 
conditions or resource states may be involved in structuring the community (McGee and Eaton, 
2015). Furthermore, in niche overlap simulated communities, niche utilization values were 
randomly shuffled for each species through 50,000 interactions. The observed mean niche overlap 
and the mean overlap in the simulated communities, the probability of the observed mean overlap 
and the standard effect size (SES) were calculated. The hypotheses stated that niche overlap does 
not differ significantly from the expected level of overlap, even in the absence of species 
interactions. In general, the species partitioned the niche (an assemblage that was structured by 
resource partitioning) when the overlap value was significantly lower than the expected level of 
overlap, even in the absence of species interactions. In general, the species shared the niche when 
niche overlap was greater than the expected level of overlap, even in the absence of species 
interactions (Miguel et al., 2016). 


3 Results 


3.1 Species diversity, biomass and soil quality changes 


In our study sites, plant coverage increased within the exclusion (Fig. 1). Species diversity 
significantly decreased (P<0.05), herbage mass significantly increased (P<0.05) in grass while 
significantly decreased (P<0.05) in Leguminosae, Compositae and weeds within the exclusion 
(Table 1). The main improvements in soil quality occurred in the topsoil (0-10 cm) within the 
fenced grassland (S1 Table). In general, grazing exclusion may be an efficient way to restore 
grassland but this practice may sacrifice biodiversity. 


Table 1 Changes in species diversity and biomass inside and outside the fence 


Vegetation property Index Outside the fence Inside the fence t-test 
Diversity index R 13.73+1.87 10.73+1.58 P=0.00 
H' 10.14+1.21 7.41+1.28 P=0.00 
D 8.04+0.94 5.50+1.22 P=0.00 
J 0.88+0.03 0.84+0.04 P=0.00 
ies Citeren Fly Grass 43.57£10.23 74.14£18.99 P=0.00 
Leguminosae 37.84+16.75 11.89+46.24 P=0.00 
Compositae 5.24+5.62 1.11+1.21 P=0.03 
Ruderal 37.79+11.86 18.99+11.98 P=0.00 
Total biomass 124.44+28.49 106.13+18.69 P=0.05 


Note: Plant diversity of grassland was assessed using the richness index (R), Shannon-Wiener index (H'), Simpson index (D), and Pielou 
index (J). Biomass measured by Grass, Leguminosae, Compositae and Ruderal, respectively. Differences in mean diversity between two 
sites were assessed with paired-sample t-tests. Mean+SD. 


3.2 Relative contributions of environmental filtering and spatial factors 


The total variation ([a+b+c]) explained by environmental and spatial factors outside the fence was 
greater than inside the fence. The separate components [a], [b], and [c] were different inside and 
outside the fence (Fig. 3). The pure effects of environmental variables [a] and spatial structure [c] 
turnover increased by 31.5% and by 25%, respectively, while the variation controlled by spatially 
structured environment variables [b] decreased by 53.5% inside the fence. Only the pure effect of 
environment variables [a] was statistically significant and explained a large percentage of the 
variation in diversity both inside and outside the fence (Fig. 3). The environmental variables were 
significantly related to plant diversity between the two sites (S2 Table). For example, the silt 
content explained 81.0% of total variation in diversity outside the fenced grassland; whereas, the 
coarse sand and silt content contributed 63.6% and 15.3% of total variation in diversity inside the 
fenced grassland, respectively. Outside the fenced grassland, the plant density was significantly 
correlated with space distance (P<0.05), even after controlling for the effect of environment 
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distance in a partial Mantel test, and the correlation coefficients were higher than those inside the 
grassland (Table 2). In contrast, the correlation inside the fence was not significant (P>0.05), and 
the plant diversity displayed higher correlation coefficients with environment distance, even after 
controlling for the space distance (Table 2). 


Table 2 Mantel test and partial Mantel test correlations for community diversity, space distance, and 
environmental dissimilarity inside and outside the fence 


. Outside the fence Inside the fence 
Matrices 
R P R IP 
Environments distance 0.0094 0.423 0.078 0.065 
Environments|space distance* 0.0095 0.394 0.061 0.161 
Space distance 0.1372 0.002 0.056 0.050 
Spacelenvironments distance? 0.1370 0.001 0.060 0.072 


Note: *, plant community dissimilarity with environments distance, controlling for spatial distance; >, plant community dissimilarity 
with spatial distance, controlling for environments distance. 


(a) Outside the fence 


[a] [a] 
0.378 0.497 


P<0.001 P<0.001 


Residuals=0.137 Residuals=0.270 


Fig. 3 Variation partition results for diversity inside and outside the fence. [a]=variation explained by the 
environmental factors independent of space, [b]=variation explained by the spatially structured environment 
factors, and [c]=variation explained by spatial factors independent of any environmental factors. 


3.3 Pattern of species co-existing 


The EcoSim 7.0 analysis (Table 3) showed that the community was non-randomly structured 
inside the fence, and there was an inter-specific competition between members of this community, 
suggesting a competitively structured community (SES=2.29). It also showed that the observed 
niche overlap was significantly less than the expected overlap inside and outside the fence (SES= 
—28.17 and SES= —56.2, respectively), suggesting that the differential use of resources or 
limitation of some environmental resource helped to structure the community. 


Table 3. Co-occurrence (C-score) and niche overlap analysis of all species inside and outside the fence 


Comparison of Monte Carlo 


C-score analysis Cops and Cexp results SES Structure 
P values of obs and exp 

Inside Pobs> Pexp=0.01566 Cobs>Cexp=174.7044/50000 229 CCS 

Outside Pons? Pexp=0.27968 Cobs> Cexp=177.4621/50000 055 NCCS 

Nich overlap Comparison of Monte Carlo Ghende precnie SES Scns 

analysis P values of obs and exp 

Inside Povs<Pexp=0.0000 Obs<exp=0.75063/50000 —28.17 PLR 

Outside Pobs<Pexp=0.0000 Obs<exp=0.75082/50000 56.20 PLR 


Note: The C-score and niche overlap values of the observed community data (obs) were compared to those of the Monte Carlo 
generated pseudo-community (exp), and the number of times out of 5000 iterations and the probability (P value) that one was greater or 
lower than other were calculated. From this a standardized effect size (SES) was determined, and whether the plant communities were 
non-competitively or competitively structured (NCCS or CCS) and were not partitioning or partitioning some limited resource (PLR or 
NPLR). 
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4 Discussion 


The deterministic processes, such as environmental filtering and competitive exclusion, are the 
primary processes in shaping plant communities in the harsh environment of semi-arid grasslands 
(Laliberté et al., 2013; Wang et al., 2013). For instance, individual interactions and environmental 
filtering jointly determine community structure at the local community scale (Laliberté et al., 
2013). As expected, the most significant contribution by an environment variable was the 
important process in community assembly inside and outside the fence (Fig. 3). Among the 
examined environment factors, soil particle size distribution, such as silt content outside the fence 
and coarse content inside the fence, explained a significant fraction of the total variance that was 
detected in the plant community assembly based on redundancy analysis (S2 Table). Grazing 
exclusion increased the fine sand and silt+caly content in the topsoil (Su et al., 2005; Wu et al., 
2014). Moreover, the fenced management improved the soil structure and soil-water capability 
(Wu et al., 2014). Soil particles should be well correlated with soil hydrology physical processes. 
For example, soils that have large particles and pores can hold large amounts of water but do not 
retain the water since the water leaches through the pores. In contrast, clay and silt soil have low 
porosity and high water retention (Wu et al., 2014). Therefore, soil particle size distribution may 
play an important role in determining the soil water availability, which largely determines the 
structure of community assembly in arid and semi-arid habitats (Chesson et al., 2004; Miranda et 
al., 2011). Our finding showed that there were more silt and clay in topsoil inside the fence in 
comparison to outside the fence (S1 Table), resulting in a higher water retention in the topsoil and 
an increase in the vertical segregation of soil water availability. In other words, variation in the 
vertical soil profile may create a spatially heterogeneous water resource. 

Niche overlap value of the observed community data was significantly less than the expected 
value inside and outside the fence (Table 3), suggesting a differential use of resources or 
limitation of some environment resource that could be in help to structure the community. 
However, both types of processes (environmental filtering and competitive exclusion) did not 
work together to determine species diversity and coexistence outside the fence (Fig. 3; Table 3). 
Disturbance has been acknowledged as one of the important factors affecting species assembly by 
creating ecological niches, changing soil processes, and leading to competition (Grime, 2006; 
Bello et al., 2013). Therefore, free grazing by sheep outside the fence may cause a dilutive effect 
of competition among different species, especially between annual species and perennials. For 
instance, selective foraging exerts a local decrease in the chance of a competitive exclusion of 
species, thereby allowing the subordinate species to coexist and increasing species richness 
(Grime and Mackey, 2002; Bello et al., 2013; Nathan et al., 2016). Grazing also introduces more 
heterogeneity into grasslands, which may be beneficial for structural complexity and diversity 
(Jerrentrup et al., 2015). Heterogeneous structures may allow different species to occupy different 
microhabitats and provide ecological niches promoting biodiversity (Mayfield et al., 2010). This 
was confirmed as the spatially structured environmental variation explained 47.7% of species 
diversity; whereas its explanatory power (22.2%) decreased inside the fence (S2 Table). In 
contrast, grazing exclusion led to an increase in relative biomass of grass, while the biomass of 
Leguminosae, Compositae and ruderal decreased. The transition of above-ground productivity 
may reflect competition among different families (Table 1), suggesting that the grass family has a 
greater competitive ability for space and resources. Therefore, competition between co-existing 
species plays a significant role inside the fence (Table 3). 

The different processes, such as niche and dispersal processes, may be not mutually exclusive 
and jointly determine community assembly across different scales, regions and vegetation types 
(Wang et al., 2013; Li et al., 2015). Contrary to this viewpoint, however, the community 
composition in this study was controlled by a negligible effect of space that was independent of 
environment (Fig. 3). Seed limitation strongly influences plant community assembly in grassland 
(Marteinsdóttir, 2014). The dispersal filter was not acting in a semi-arid grassland system 
dominated by a tussock-forming Stipa tenacissima, because there was a dense and evenly 
distributed soil seed bank (Luzuriaga, 2012). In our study region, mean annual wind velocity was 
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2.8 m/s and the prevailing winds mainly occurred in April and May. If we assume that the seed 
bank was evenly distributed in this region, the climate condition may have been beneficial for 
seed dispersal. In addition, patches of bare soil created by tramping may have offered space for 
propagation and dispersal of seeds (Luzuriaga, 2012). 


5 Conclusions 


Our findings indicate that in semi-arid grassland ecosystems dominated by tussock-forming S. 
breviflora (perennial grass), rapid changes in topsoil quality that occur following grazing 
exclusion were associated with a fast recovery of biomass production in grass. The different 
processes, environmental filtering and competitive exclusion processes, were not mutually 
exclusive and collectively determined community assembly after grazing the exclusion. Shifts in 
the relative importance of different processes may be due to certain ecological changes. First, 
grazing exclusion increased the fine sand and silt+clay content in the topsoil, which may correlate 
with soil hydrology and influence the soil water availability. Second, a rapid increase in the 
relative biomass of the grass family may enhance its competitive ability for space and resources. 
It can be inferred that grazing exclusion may be an efficient way to restore grassland but this 
practice may sacrifice biodiversity. For the maintenance and restoration of diversity, designers 
should be aware of this risk and account for the steering community assembly to maximize 
species diversity, especially when manipulating processes that introduce new species and promote 
environmental heterogeneity. 
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Appendix 
S1 Table Soil organic carbon content and particle size distribution to a depth of 40 cm inside and outside the 
fence 
Soil depth 
Parameter Sites 
0-5 cm 5-10 cm 10-20 cm 20-40 cm 
Outside the fence 0.27+0.07 0.27+0.09 0.4+0.12 0.32+0.12 
Soil organic carbon (%) Inside the fence 0.29+0.07 0.3440.22 0.37+40.11 0.27+0.08 
t-test (P value) 0.41 0.27 0.12 
Outside the fence 2.06+0.48 2.9240.73 4.96+1.35 5.36+0.92 
Clay (<0.002 mm) (%) Inside the fence 2.44+0.69 4.36+1.85 23 ES 2.34+0.96 
t-test (P value) 0.10 0.01 0.00 
Outside the fence 18.52+3.63 25.3746.75 41.41+9.82 44.22+8.09 
Silt (0.002—0.05 mm) (%) Inside the fence 19.7944.88 30.9547.99 44.31+4.09 19.8347.4 
t-test (P value) 0.46 0.06 0.00 
Outside the fence 35.3541.21 31.75£2.33 24.94+4.11 23.743.47 
Fine sand (0.05—0.1 mm) (%) Inside the fence 32.4343.2 28.16+3.74 22.94+2.18 34,543.63 
t-test (P value) 0.01 0.00 0.00 
Outside the fence 44.06+43.5 3997509 28.69+6.96 26.72+5.72 
Coarse sand (0.1-2 mm) (%) Inside the fence 45.34+4.45 36.53+6.03 25.51+3.33 43.33+5.77 
t-test (P value) 0.38 0.15 0.00 


Note: The difference between soil organic carbon (SOC) contents wasn't statistically significant inside versus outside the fence. SOC 
inside the fence increased at depths of 0-5 and 5-10 cm, and declined at depths of 10-20 and 20-40 cm when compared with outside 
the fence. The soil particle contents inside the fence, such as clay (<0.002 mm) content, was significantly (P<0.05) increased at depth of 
0-20 cm but decreased at a depth of 20—40 cm. The silt content at depths of 0-5, 5-10 and 20—40 cm exhibited differences similar to 
those of the clay contents inside the fence. The fine sand content decreased significantly (P<0.05) at depths of 0-5 and 5-10 cm but 
increased significantly (P<0.05) at depth of 20-40 cm compared to outside the fence. Inside the fence had higher coarse sand content 


(0.1-2 mm) at all depths than the outside, but only at depth of 20—40 cm was statistically significant (P<0.05). 


S2 Table Explanatory environmental and spatial variables selected by the forward selective procedure in the 
redundancy analyses inside and outside the fence 


Sites Variables Order Variables R? R’Cum Adj R?Cum F P value 
1 Silt 0.810 0.810 0.808 417.013 0.001 

Environment 2 Coarse sand 0.006 0.852 0.849 28.039 0.001 

3 Cover 0.043 0.859 0.854 4.417 0.035 

1 PCNM4 0.137 0.137 0.129 15.614 0.001 

2 PCNM20 0.069 0.206 0.190 8.389 0.003 

3 PCNM9 0.061 0.267 0.244 7.950 0.005 

OE 4 PCNM8 0.059 0.325 0.297 8.239 0.007 
the fence 5 PCNM2 0.037 0.362 0.329 5.481 0.022 
Space 6 PCNM10 0.036 0.399 0.360 5.639 0.019 

7 PCNM21 0.033 0.432 0.388 5.278 0.015 

8 PCNM17 0.029 0.460 0.413 4.868 0.018 

9 PCNM28 0.029 0.489 0.438 5.036 0.024 

10 PCNM40 0.028 0.517 0.463 5.241 0.016 

11 PCNM14 0.025 0.542 0.485 4.722 0.030 

1 Coarse sand 0.636 0.636 0.633 171.536 0.001 

Inside i 2 Silt 0.153 0.789 0.785 70.099 0.001 
the fence PaviTonment 3 Cover 0.010 0.799 0.793 4.782 0.022 
4 Fine sand 0.012 0.811 0.803 6.138 0.015 


To be continued 
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Continued 
Sites Variables Order Variables R R?Cum Adj R?Cum F P value 
1 PCNM25 0.070 0.070 0.061 7.426 0.006 
2 PCNM35 0.062 0.133 0.115 6.952 0.003 
Linse Space 3 PCNM39 0.051 0.183 0.158 5.973 0.015 
the fence 
4 PCNM13 0.046 0.230 0.197 5732 0.015 
5 PCNM6 0.041 0.271 0.233 5.342 0.020 


Note: We used forward selection for all environmental and PCNM variables to obtain the significant environmental and PCNM 


variables, respectively. Outside the fence, the silt, coarse sand content, and the cover of S. breviflora were the main environmental 
variables; and eleven PCNM variables (#2, 4, 8, 9, 10, 14, 17, 20, 21, 28, and 40) were selected as the main spatial variables. Together, 
these explanatory environmental and PCNM variables explained 85.4% and 48.5% of the total plant diversity variances, respectively. 
Inside the fence, the coarse sand, silt, fine sand content, and the cover of S. breviflora were the major environmental factors extracted; 
and five PCNM variables (#6, 13, 25, 28, 35, and 39), representing the major spatial predictor were generated. These explanatory 


environmental and PCNM variables accounted for a total of 80.3% and 23.3% of the plant diversity variance, respectively. 


